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1 Introduction 

The origin of sub-eV neutrino masses is a central issue in particle physics. As right-handed 
neutrinos do not exist in the Standard Model (SM), which thus predicts massless neutrinos, 
new particles are necessary to explain neutrino masses [1], e.g., gauge singlet fermions 
in the Type I [2-7] Seesaw Mechanism, or scalar and fermionic SU(2)i triplets in the 
Types H [8-11] and HI [12] scenarios. Searches for these degrees of freedom constitute 
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Figure 1. (a) Leading order Feynman diagram for pp —?> T+T production, (b-h) 

0{as) corrections. Drawn with JaxoDraw [32]. 


an important component of hadron collider programs; see Refs. [13-16] and references 
therein. Furthermore, the maturity of the formalism underlying QCD corrections in hadron 
collisions, which are required for predicting accurate production rates and distribution 
shapes, readily permit their application to beyond the SM (BSM) processes. 

For heavy Seesaw partners with sub-TeV masses, the dominant hadron collider produc¬ 
tion mode is through the Drell-Yan (DY) charged current (CC) and neutral current (NC) 
processes [17-23], shown in figure 1(a). For TeV-scale systems and above, the fusion 
channel becomes dominant [19, 24, 25]. A catalog of resonant Seesaw partner production 
modes in hadron collisions is given in Ref. [19]. 

In the Type I Seesaw, production cross section is known at next-to-leading order (NLO) 
in QCD [26] and estimated at next-to-next-to-leading order (NNLO) via a iF-factor* [25]. 
For the Type II case, rates are known at NLO [27, 28], NLO with next-to-leading logarithm 
(NLL) recoil and threshold resummations [29], and automated at NNLO [30, 31]. 

Pair production of heavy Type III Seesaw leptons has, until now, been evaluated only 
to leading order (LO) accuracy. For rriT = 100 GeV —2 TeV, we report the NLO iL-factors: 


1.17-1.37 

at 


13 TeV, 

(1.1) 

1.17-1.36 

at 

V^ = 

14 TeV, 

(1.2) 

1.14-1.29 

at 

V^ = 

100 TeV, 

(1.3) 


with scale uncertainty of at 14 TeV and ±1% at 100 TeV, and are comparable to other 
DY-type processes in Seesaw models. The NLO differential 7L-factors^ for heavy lepton 
kinematics are largely flat for TeV-scale mp, suggesting that naive scaling by the total 
rNLO jg reasonably justified. 

In this study, production rates of TeV-scale Type III Seesaw lepton pairs at 0{as) accu¬ 
racy are presented for pp collisions at ^/s = 13, 14, and 100 TeV. Differential distributions 

‘The N^LO A'-factor is defined as K = , where is the lowest order (m = 0), or Born, 

cross section and is the N’"LO-corrected cross section. 

^The differential N^LO Y-factor with respect to observable O is defined 

(da^'"^°/dd'j / {da^°/ddy 
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as 














at NLO and NLO with leading logarithm (LL) resummation of TeV-scale lepton kinematics 
are presented for the first time at 14 TeV. The fixed order (FO) calculation is carried out 
via phase space slicing (PSS) [33-36]. The calculation of the dilepton system’s transverse 
momentum, qt, follows the Collins-Soper-Sterman (CSS) formalism [37-39]. This text con¬ 
tinues in the following order: In section 2, we summarize the Type III Seesaw model and 
comment on experimental constraints. The PSS and CSS formalisms are briefly introduced 
in section 3. Due to the factorization properties of DY-type systems, the fixed order and 
resummed results reduce to convolutions over tree-level quantities; technical details are 
relegated to appendices A and B. Results are reported in section 4. We summarize and 
conclude in section 5. 

2 Type III Seesaw Mechanism 


2.1 Model Lagrangian 

The Type III Seesaw [12] generates tree-level neutrino masses via couplings to SU(2)2, 
triplet leptons with zero hypercharge. In terms of Pauli matrices a°‘, the left-handed (LH) 
fields are denoted by 


Sl = 


( Si ^ Sj T i^l 

\V2^l -Siy’ ^ V2 


( 2 . 1 ) 


where have U(1)em charges Q = ±1, and the right-handed (RH) conjugate fields are 


-Si’ J ' 


( 2 . 2 ) 


Chiral conjugates are related by = {iP^)r = {'•PlY, where Yl/r = Pl/rY = ^(1 T 1^)Y- 
For a single generation (but generalizable to more), the model’s Lagrangian is 


71 Type III = TIsm -\- Ct + Py, 


(2.3) 


where UsM is the SM Lagrangian, the triplet’s covariant derivative and mass are given by 
£t = ^ TV IY^l] - (^^Si^ + mT^S+‘’ + H.c) , (2.4) 

and the SM LH lepton (L) and Higgs ($) doublet fields couple to S^ via the Yukawa 
coupling 

Cy = -vtL Si + H.C., v = V2(4>) « 246 GeV. (2.5) 

Dirac masses are then spontaneously generated after electroweak symmetry breaking (EWSB): 

Cy ™ Cy = -^{v + h)l7ZCY - yriv + h)eZC+^ + H.c., (2.6) 
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leading to the neutral fermion mass matrix 




S + H,c^ = - 


—f 0 yTv/V^\ f V' 




.yTv/y/2 rriT J 



+ H.C. (2.7) 


The Seesaw mechanism proceeds by supposing rriT 2/r(^)) leading to light/heavy mass 
eigenvalues 

„ 2,,2 

and mheavy ~ mr- (2.8) 


flight 


yW 

2mT 


Thus, tiny neutrino masses follow from mixing with heavy states, whereby light (heavy) 
mass eigenstates align with the doublet (triplet) gauge states. For yT comparable to the 
electron’s SM Yukawa, sub-eV mught can be explained by sub-TeV niT, a scale within the 
LHC’s kinematic reach. 

We combine the S helds and their conjugates into physical Dirac and Majorana helds: 


r- = S" + f + = f-^, f ° = + sic (2.9) 

In the gauge basis and in terms of T, the triplet interaction Lagrangian Ct is written as 
^Gauge Basis ^ f- ^ 1 fO 

- {eA^r + gcos0wZ^,j>^) T" - - gT0W+^>^f-(2.10) 

Our aim is to report the 0{as) corrections to heavy lepton pair production, which are 
independent of the mixing between the gauge states T and mass states T, i [i = 
and (w- = 1)2,3). For the remainder of the study, we generically denote the mixing as 
Y and e, and write 


= Y +£ i^, r° = yr° + ei/m, |Y|~0(1), |e| < l. (2.11) 


The resulting interaction Lagrangian in the mass eigenbasis relevant to our study is 
^Mass Basis ^ _ ^t) 

- ^ (eYA^^^ + g cos QwYZ^^^) T" - gY'^W-^^T^ - gY*'^W^^^(^A2) 

2.2 Constraints on Type III Seesaw Lepton Production 

For a review of constraints and phenomenology of the Type III Seesaw, see Refs. [15, 40-45]. 

• Collider Production and Decay: CMS experiment searches for production 

and decay into t^l ~have restricted the cross section and branching ratio to [46] 

a{pp r°r^) X BR(r°r± ^ i+re'^ ^t) < 12 fb at the 95 % c.l. ( 2 . 13 ) 


For equal doublet-triplet mixing among the SM leptons, this translates to the bound 


ttit < 278 GeV at the 95% C.L. (2-14) 

With 20.3 fb“^ of 8 TeV LHC data, searches carried out by the ATLAS experiment for 
—>■ 4jT^ excludes [47], depending on mixing parameters, 

rriT < 325 — 540 GeV at the 95% C.L. (2.15) 
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Throughout this study, we take and to be mass degenerate. Electroweak (EW) 
corrections at one loop induce a mass splitting of Auit ~ 160 MeV for rriT > 100 GeV [42, 
48, 49], and is thus negligible. Eor differential distributions, we use representative mass 

ruT = 500 GeV. (2.16) 

As in the LO case, the total partonic and hadronic cross sections at NLO and NLO+LL 
factorize into a product of the mixing parameter |y| and a mixing-independent “bare” cross 
section uq: 

^ |y|2 ^ a^LO^p (2.17) 

Therefore, we express our results in terms of cjo and do not choose any particular |y|. 
Furthermore, factorization implies that that total and differential NLO iL-factors are inde¬ 
pendent of |y|. 

3 Heavy Lepton Pair Production at 0{as) in Hadron Collisions 

Here we outline the PSS [33-36] and GSS [37-39] formalisms, which we use to calculate the 
processes 


pp^W^*^T°T^ and p p-f* / Z* ^ T+T-, (3.1) 

at NLO in QGD and the transverse momentum qt of the dilepton systems at LL. With 
PSS and CSS, the inclusive NLO and NLO+LL results factorize and can be expressed in 
terms of tree-level, partonic cross sections. Such technical details are given in appendices 
A and B. For simplicity, we generically denote processes in Eq. (3.1) and their radiative 
corrections by 


p p ^ T T and p p T j. (3.2) 

We note that these corrections are not unique but are well-known and general for the 
production of any SU{2)l triplet color-singlet, e.g., [50, 61[. However, unlike previous 
studies, we investigate the 0{as) effects on the kinematic distributions of TeV-scale leptons. 

3.1 Phase Space Slicing 

To evaluate TT production at NLO, we follow the usual procedure: evaluate virtual and 
radiative corrections to the LO process in d = 4 — 2e dimensions; collect soft divergences, 
which cancel exactly; collect collinear divergences, which cancel partially; and subtract 
residual collinear poles from parton distribution functions (PDFs). 

For an n-body LO process, we divide, or slice, the phase space of its (re + l)-body 
correction into soft and collinear kinematic regions. For radiation energy Ej, partonic 
c.m. energy \/I, and small dimensionless cutoff parameters 5s, 5c, a volume of the (re+ 1)- 
body phase space is soft if 

Ej < ^5s. (3.3) 
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For partonic-level invariant masses and momentum transfers 

Sik = (Pi + Pkf and iik = [pi - Pkf, (3.4) 

where indices i, k run over initial- and final-state momenta, a region of phase space is 
collinear if 

Sik, \iik\ < s X 6c- (3.5) 

A volume is hard {non-collinear) if not soft (collinear). Exact choices of 6s, 6c do not 
matter: dependences on 6s, 6c cancel for sufficiently inclusive processes [36]. However, so 
soft and collinear factorization remain justified, one needs 


(fc < <^5 < 1. 

The hard-non-collinear TTj process is then finite everywhere and given by 


(3.6) 


(^i3){PP^T T j X) = 

Y. [ [ ^^2 [fa/p{Ci,pj)fb/p{^2, p)) + (1 o 2)] a^{ab ^TT j), (3.7) 


a,b=q,q',g 


where is the Born-level TTj partonic cross section. For a,b £ {<l,Q',g} with q £ 
{u,d,c,.s}, the PDF fa/p{^i, p'j) is the likelihood of parton a carrying away longitudinal 
momentum fraction from proton p evolved to a factorization scale pj. The c.m. beam 
energy ^/s and partonic c.m. energy are related by s = and we denote the threshold 

at which TT production occurs by tq\ 


s {ruT + mip)^ 

tq = mm - = - 

s s 


(3.8) 


In the soft/collinear limits, amplitudes for soft, soft-collinear, and hard-collinear radi¬ 
ation factorize into divergent expressions proportional to the (color-connected) Born am¬ 
plitude. The poles are grouped with virtual corrections and the PDFs, resulting in a finite 
expression given by [36] 


o'(2)(pp-t T T X) = 

fjHC Y [ [/a/p(6,/^/)/6/p(6,/^/) + (1 ^ 2)] (T(2)(a6^ T r)(3.9) 

a,b=q,qdT0 ^ro/Ci 

fT(2) — (7 +(J (7 +(J , (o.lUj 

were d® is the Born-level TT partonic cross section, is its 0{as) virtual correction, 
and are the soft and soft-collinear radiation terms, and is the hard-collinear 
radiation correction. Inclusive triplet lepton production at NLO is now reduced to a sum 
of two- and three-body processes: 

^TTX) = a^2){pp ^TTX)+ a^s){pp ^TTj X). (3.11) 
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Figure 2. The 14 TeV LHC (a) T'^T^ and (b) r+T“ production cross section, divided by |Fp, at 
LO (solid) and NLO (dash) as a function of 5s with the two- (dot-upside down triangle) and three- 
body (dot-triangle) NLO contributions. Panel: Ratio of NLO to LO with Monte Carlo uncertainty. 


Using the inputs from section 4 and representative triplet mass rriT = 500 GeV, in 
figures 2(a) and 2(b) we show as a function of soft cutoff 5s the 14 TeV and T~^T~ 

bare production cross sections at LO (solid) and NLO (dash) [Eqs. (3.11)] with the two- 
body (dot-upside down triangle) [Eqs. (3.9)] and three-body (dot-triangle) [Eqs. (3.7)] NLO 
contributions. In the panel, we show the NLO iL-factor with Monte Carlo uncertainty. The 
negative value of cr( 2 ) for 5s < 0(0.5) is due to large hard-collinear PDF subtractions. The 
NLO result is insensitive to 5s for 5s < 1 — 3 x 10“^, reflecting the large but fine cancellation 
of Eqs. (3.7) and (3.9). For example: for = 1 x 10“^, the three- and two-body calculations 
are approximately -1-10.2 and —9.0 times the LO cross section. For ^5 > 3 x 10“^, the K- 
factor plummets, indicating the importance of terms proportional to powers of 5s in the 
two-body expression, and hence a breakdown of soft/collinear factorization. 


3.2 Collins-Soper-Sterman Transverse Momentum Resummation 

As a color-singlet process, colored initial state radiation (ISR) is the dominant contribution 
to the TT system’s qt spectrum. For dilepton invariant mass My*, the distribution has 
the power series 


dajpp ^TT X) 

dq^ 


E 

k=l 


Ak a^(M2*) log(2^-i) 


Qt 


Ak ~ 0(a^ 


(3.12) 


and indicates a breakdown of the perturbative description in the (^/My-t —)• 0 limit. The 
leading a* log(My*/g|.) logarithm in FO calculations is only reliable when [39] 


log 


My* 


A, 


QCD 


;:$> 1 and log' 


My* 

qr 


^ log 


My* 


A, 


A, 


QCD 


QCD 


= 0.2 GeV. 


(3.13) 
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That is, as{Myt,) must be perturbative and the scales associated with the process must 
be comparable. As the Seesaw mass scale is pushed higher [46, 47], so too does the scale 
at which artihcially large logarithms appear. For rriT = 500 GeV (1 TeV), FO predictions 
breakdown at qx ~ 55 (95) GeV, and resummation of recoil logarithms become necessary 
to describe qx below this threshold. 

Fortunately, as qx/My* —)• 0, gluon radiation factorizes. This permits one to reorganize, 
sum, and exponentiate large logarithms in Eq. (3.12), resulting in an all-orders expression 
in terms of the Born process. The resummed distribution with respect to q^, Tfy*, and 
dilepton rapidity y, is [39] 

The integral is over the impact parameter b and is the Fourier transform of qx] the zeroth 
order Bessel function Jq emerges as a simplification. W expands to a perturbative and non- 
perturbative set of universal Sudakov form factors, and a process-dependent luminosity 
weight W: 

W = W. (3.15) 

Expressions for Sjyp and Sp are in appendix B, in Eqs. (B.IO) and (B.14). For partons i, k, 
W is 


w= 

i,k=q,q',g 


+ (6 ^ 6 ) 


(3.16) 


are the Fourier transformed transverse momentum-dependent (TMD) PDFs evolved to 
impact scale b and collinear factorization scale fif. To LL, W is given in Eq. (B.25). 

The resummed result describes well the qx My* behavior due to the Sudakov sup¬ 
pression. But because of neglected terms proportional to powers of qx, it underestimates the 
spectrum at qx ^ My*, precisely where the FO calculation becomes reliable. To describe 
accurately qx everywhere, one introduces the auxiliary function that matches the 

asymptotic FO (resummed) behavior at small (large) qx- Gombining the three expressions, 
the total, matched spectrum is given by [51, 52] 

daMatched ^ ^^Resum ^ ^fo-Asymp 

dqx dqx dqx dqx 

The area bound by the c^^jAlatched jg then normalized to the total rate of 

Eq. (3.11) [67]. Individual terms of Eq. (3.17) are given in Eqs. (B.6), (B.26), and (B.27). 


4 Results 

Tree-level results are calculated using helicity amplitudes. The Cuba library [53] is used for 
Monte Carlo integration; numerical uncertainty is negligibly small. Events are output in 
Les Houches Event (LHE) format [54]. Rates and shapes are checked by implementing the 
Lagrangian of Eq. (2.12) into FeynRules 2.0.6 [55, 56] and using MadGraph_aMC@NLO 









v5.2.1.0 [57] (MGS). Rates are also in agreement with literature [42]. We take as SM 
inputs [58] 

Mz = 91.1876 GeV, = 1/127.944, sin|jg(0vi/) = 0.23116. (4.1) 

The Cabbibo-Kobayashi-Masakawa (GKM) matrix is taken to be unity, introducing a 
percent-level error that is no larger than the estimated O(a^) contributions. The CJ12mid 
NLO parton distribution functions (PDFs) [59] with as{Mz) = 0.118 are used. Since 
2m'];' > M\y,Mz, the factorization (fif) and renormalization (/r^) scales are fixed to the 
sum of the triplet lepton masses 


^0 = fif = = 'niT + ~ ( 4 - 2 ) 

as{Mz) is run to /.io at one-loop in QGD with nf = 5. Setting Uf = 6 increases the three- 
body channel by 2 — 3%, but the total NLO cross section by less than For total cross 

section calculations, we choose soft and collinear cutoffs 

(55 = 1.0 X 10"^ and dc = (55/100. (4.3) 

PSS involves fine cancellation of large numbers; for differential distributions, 6s is relaxed 
to 

6s = 3.0 X 10"^. (4.4) 

Born-level TTjX events can be generated efficiently by implementing soft and collinear 
cuts Eqs. (3.3) and (3.5) into MG5; see appendix A.3. Similarly, TTX events (sans the 
hard-collinear PDF subtraction) can be efficiently produced by applying an appropriate 
scaling. 

For plots in this section, the LO (NLO) curve is denoted by a solid (dashed) line. 

4.1 Total T°r=*= and T+T” Production at NLO 

Figures 3(a) and 3(b) show, respectively, the total CG and NC production cross section, 
divided by the mixing parameter \Y\‘^, as a function of heavy lepton mass for ^/s = 14 and 
100 TeV. For ttt-t = 100 GeV — 2 TeV, the NLO production rates range 


7 ab — 24 pb at 14 TeV (4-5) 

10 fb - 230 pb at 100 TeV. (4.6) 

The corresponding NLO T^T~ cross sections range 

5 ab — 12 pb at 14 TeV (4-7) 

5 fb - 130 pb at 100 TeV. (4.8) 


The CG rate is approximately twice as large as the NC channel due to W boson charge 
multiplicity. In the low-(high-)mass range, transitioning from 14 to 100 TeV increases the 
total cross section roughly by a factor of 10 (1000) for both processes. 
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Figure 3. The 14 and 100 TeV total (a) and (b) T+T production cross section, divided 

by |Tp, at LO (solid) and NLO (dash) as a function of rriT. Panel: Ratio of NLO to LO. 

The panels of figure 3 show the NLO iL-factor at 14 and 100 TeV. For T^T^, the ratios 
span 


1.17-1.36 at 14 TeV (4.9) 

1.14 - 1.17 at 100 TeV. (4.10) 

For T^T~, they range 

1.19 - 1.33 at 14 TeV (4.11) 

1.15 - 1.29 at 100 TeV. (4.12) 

At lower collider energies, AT-factors are larger for heavier niT due to the rarity of anti¬ 
quarks possessing sufficiently large momentum at LO. At NLO, this is compensated by 
large Bjorken-x gluons undergoing high-py’ 9 ^ Q splitting. CC and NO iL-factors are 
appreciable and, due to their color structures, comparable to those of the Seesaw Types 
I [25] and II [27-31]. Table 1 summarizes these results for representative mj' at y/s = 13, 
14, and 100 TeV pp collider conhgurations. 

4.2 Scale Dependence 

Higher order QCD corrections are necessary to further reduce theoretical uncertainty. To 
quantify and estimate the size of these contributions, the default scale po = 2mT is varied 
over the range 

po X 0.5 < R < Ro X 2. (4.13) 

The scale variation is then defined as the ratio of the NLO rate evaluated at scale p to the 
same rate at p = pQ. Figures 4(a) and 4(b) show, respectively, the scale variation band of 
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Table 1. TT NLO Cross Sections and NLO if-Factors for Various pp Collider Configurations 




13 TeV 

14 TeV 

100 TeV 

rriT [GeV] 

^7VLO/|y|2 

j^NLO 

^NLO/\y\2 [fb] 

j^NLO 

^AfLO/|y|2 [fb] 

j^NLO 

100 

21.7 X 10^ 

1.17 

23.9 X 103 

1.17 

229 X 103 

1.17 

300 

370 

1.17 

423 

1.17 

7.07 X 103 

1.15 

500 

43.7 

1.19 

51.9 

1.19 

1.34 X 103 

1.14 

700 

8.75 

1.21 

10.8 

1.20 

433 

1.14 

900 

2.22 

1.23 

2.89 

1.23 

182 

1.14 

1000 

1.18 

1.25 

1.58 

1.24 

125 

1.14 

1500 

62.4 X 10-3 

1.31 

97.8 X 10-3 

1.30 

28.8 

1.14 

2000 

3.60 X 10-3 

1.37 

6.96 X 10-3 

1.36 

9.58 

1.14 

rp+J^- 


13 TeV 

14 TeV 

100 TeV 

rriT [GeV] 

^7VLO/|y|2 [fb] 

j^NLO 

^iVLO/|y|2 [fb] 

j^NLO 

^AfLO/|y|2 [fb] 

j^NLO 

100 

11.2 X 103 

1.20 

12.4 X 103 

1.20 

127 X 103 

1.29 

300 

184 

1.18 

211 

1.18 

3.70 X 103 

1.20 

500 

21.4 

1.19 

25.5 

1.19 

692 

1.17 

700 

4.26 

1.21 

5.30 

1.20 

221 

1.16 

900 

1.10 

1.23 

1.42 

1.22 

92.3 

1.15 

1000 

588 X 10-3 

1.24 

782 X 10-3 

1.23 

63.5 

1.15 

1500 

35.7 X 10-3 

1.29 

53.7 X 10-3 

1.28 

14.4 

1.15 

2000 

2.74 X 10-3 

1.35 

4.81 X 10-3 

1.33 

4.78 

1.15 


the CC and NC processes at NLO as a function of rriT', the lower panels show the NLO 
iC-factor for the three scale choices. The default scale is denoted by a solid line at 1. The 
high (low) scale scheme is denoted by right-side up (upside down) triangles and is found 
to decrease (increase) the total cross section. This suggests that the renormalization scale 
evolves to smaller values faster than the factorization scale evolves PDFs to larger 

values, and also leads (accidentally) to a vanishing scale dependence for ttt-t ~ 100 GeV. 
The behavior is consistent with other TeV-scale Seesaws mechanisms [25]. For the rriT 
studied, the CC and NC calculations exhibit maximally a scale dependence at 14 TeV; 
this reduces to ±1% at 100 TeV for the same mass range considered. The 14 TeV NLO 
iL-factor varies maximally for the CC process and for the NC process. The slightly 
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(a) 


(b) 


Figure 4. Scale dependence of (a) and (b) T+T production at NLO as a function rriT- 


Table 2. 14 TeV Scale Dependence of TT Production Cross Sections at NLO 


Channel 

Scale Choice p 

mx [GeV] 

500 

1000 

1500 

2000 

2000 (100 TeV) 

'J'O'J'i: 

0.5 X pq 

+2.3% 

+3.0% 

+4.1% 

+5.3% 

+1.3% 

2x po 

-2.1% 

-3.5% 

-4.5% 

-5.6% 

-1.3% 

X+rp- 

0.5 X pq 

+2.0% 

+3.1% 

+3.7% 

+4.7% 

+1.4% 

2x Po 

-2.3% 

-3.3% 

-4.4% 

-5.6% 

-1.2% 


larger scale dependence in the itT-factors than the total cross sections is due to the larger 
scale dependence of the LO result. The size of the //-dependence suggests that 0(0^) 
effects are small, consistent with Refs. [25, 30, 31]. Results are summarized in table 2 for 
representative ttit- 

4.3 14 TeV Kinematic Distributions at NLO and NLO+LL 

Figure 5(a) shows the 14 TeV NLO differential distribution, divided by the mixing parameter 
|Tp, with respect to the px of in production. The panel shows the differential 

NLO iC-factors. At low (high) px the change at NLO is small (large) and follows from 
the T^T^j channel at 0{as). The transverse recoil of dilepton system from hard ISR 
propagates to individual leptons, thereby providing an additional transverse boost. As the 
jet energy softens or is radiated more collinearly to its progenitor, its px vanishes, and 
kinematics at NLO approach those at LO. 
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Figure 5. The 14 TeV differential distribution, divided by |Fp, with respect to (a) px and (b) y 
of in production at LO (solid) and NLO (dash). Panel: Ratio of NLO to LO. 


In figure 5(b), the rapidity (y) distribution of in production is presented. 

We observe that QCD corrections have a small impact on the rapidity distribution shape 
as indicated by the mostly flat it'-factor, with only a slight upwards bump at small y, 
where pp ^ Pz- Similar to the pp spectrum, new kinematic channels at 0(as) all involve 
high-pr ISR and do not induce longitudinal boosts, leaving the y distribution shape largely 
unchanged. 

Similar pp and y behavior are observed for in and in the NC processes. 

Aside from the Majorana nature of light {i') and heavy (T^) neutral leptons, and the 
relative CC and NC production rates, i.e., (T{pp —)• T^T^)/a{pp —)• T^T~) ~ 2, observing 
vector-like coupling of heavy leptons to electroweak gauge bosons is a critical test of the 
Type III Seesaw mechanism [42]. As in the SM, this done by measuring the polar distribu¬ 
tion made by, for example, in production in the dilepton rest frame with respect 

to dilepton system’s direction of propagation in the lab frame. Symbolically, the observable 
is given by 


cos 6* 


P*p-q 

\pIM' 


(4.14) 


where is the 3-momentum of lepton T in the TT frame and q = pp + p^ in the lab 
frame. Figures 6(a) and 6(b) show, respectively, the CC and NC cosO* distribution. At LO 
and NLO, the vector coupling structure is clear. The uniform iL-factor follows from 0{as) 
corrections involving only initial-state partons and amount simply to a boost of the dilepton 
system in the lab frame. The 0(as) effects are unraveled in constructing Eq. (4.14) and 
subsequently affect only the normalization. 

As discussed in section 3.2, the transverse momentum distribution qp of the 
dilepton system is ill-defined at qp My* in FO perturbation theory and requires resum- 
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Figure 6. Same as figure 5 but with respect to cos0*, as defined in Eq. (4.14), of (a) in 
production and (b) T~^ in T+T production. 


mation of logarithms associated with soft gluon radiation. For representative masses (a) 
rriT = 500 GeV and (b) 1 TeV, figure 7 shows the various contributions to qt spectrum: 
the asymptotic (dash) and FO (solid) terms at 0{as), the resummed rate at LL (dot), and 
combination of the pieces. The lower panel shows the ratio of the combined result to the 
FO result. For rriT = 500 GeV (1 TeV), the FO calculation overestimates the combined 
differential rate for qx < 25 GeV. At qx ~ 55 (95) GeV, the estimate from Eq. (3.13), 
the FO result remains about 35% (70%) below the combined result. The largeness of the 
resummation corrections is consistent with other recoil resummations at high-scales [67]. 
The combined distributions peak at ~ 5 — 6 GeV, below which the Sudakov suppression 
from multiple soft gluon emissions overtakes the divergent nature of soft emissions. 


4.4 Discovery Potential at 14 TeV High-Luminosity LHC and 100 TeV 

In the high luminosity LHG (HL-LHG) scenario [60], detector experiments aim to each 
collect 1—3 ab“^ of data. We briefly address the maximum sensitivity to heavy triplet lepton 
pair production in this scenario. TeV-scale triple leptons decay dominantly to longitudinally 
polarized weak bosons and the Higgs [42, 44], a consequence of the Goldstone Equivalence 
Theorem, implying 

BR (r± ^ W^i^i) « 2 BR (T± ^ Zi^) « 2 BR (T^ ^ M^) « 50% (4.15) 

BR (T° ^ W+r + W-i+) « 2 BR (r° ^ + Zue) « 2 BR (T° ^ + hVi) « 50%. 

(4.16) 
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Figure 7. The asymptotic (dash) and FO (solid) at 0{as), resummed at LL (dot), and com¬ 
bined NLO+LL (dash-dot) contributions to the system’s transverse momentum differential 

distribution, qr, for mass (a) rriT = 500 GeV and (b) 1 TeV. 

For visible decay modes Z —)■ jj and h —> bb, gg, heavy lepton pairs can decay into fully 
reconstructible final-states with four jets and two high-p-p leptons that scale likep^ ~ m'r/2: 


r°T^ ^ + WZjWh W + 4j / 2j + 2b (4.17) 

T+T- ii' + ZZjZhlhh U' + 4j / 2j + 26 / 46. (4.18) 

The corresponding branching fractions are 

BR (T0r± ^ ii' + 4j/2j + 26) Ri 11.5%, (4.19) 

BR(T+T- ^ ii' + 4j72j + 26/46) « 11.6%. (4.20) 

Taking |Tp = 1, rriT = 1-5 TeV, and acceptance efficiency of .A = 0.75 [42], then after 
3 ab“^ one expects tens of heavy lepton pairs across both channels 

Ntot± = Cx cj^^O(T0r=^) X BR X ^ Ri 25.2, (4.21) 

Nt+j,- = Cx a^^'^iT+T-) x BR x ^ « 14.0. (4.22) 


To a good approximation, the kinematics of TeV -scale TT decays render the SM background 
negligible [42, 44]. Using a Gaussian estimator, the statistical significances are at the 3 —5 cj 
level: 

ct2 ’0 j'±(3 ab ^) ~ A%oj'± = 5.0, (4.23) 

ctj'+T’—( 3 ab ^) Ri ■\/Nj'+'j'— = 3.7. (4.24) 

Summing in quadrature, the combined signihcance surpasses the 6ct level (over the null 
hypothesis): 

[3 ab ^] = (Tj^0'p±[3 ab ^] 0 [3 ab ^] = 6.3. (4.25) 
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Figure 8. Required luminosity as a function of uit at 14 and 100 TeV for 5a discovery (dash-star) 
and 2a sensitivity (dash-diamond) of (a) and (b) T~^T~ production, assuming branching 

fractions and acceptance given in section 4.4. 


For rriT = 1.6 TeV, the combined significance is approximately 4.8cr, demonstrating a max¬ 
imum HL-LHC discovery potential to Type III Seesaw leptons in the rriT = 1-5 — 1.6 TeV 
range. 

Fixing the branching fractions and acceptance rates, we plot in figure 8 the required 
luminosity as a function of rriT for a 5 cj discovery (dash-star) and 2cr sensitivity (dash- 
diamond) of the (a) and (b) r+T“ channels at 14 and 100 TeV. At 14 TeV and 

after 300 (3000) fb“^, one finds sensitivity to triplet pairs up to rriT = 1-3 — 1.4 TeV (1.7 — 
1.8 TeV) in the individual CC and NC channels. At 100 TeV, we observe that with 10 fb“^ 
a 5a discovery can be achieved in the CC channel for rriT ~ 1-6 TeV and in the NC channel 
for niT ~ 1.4 TeV. At large ttit, however, taking A = 0.75 is a not justified as a different 
search methodology is required to account for the boosted kinematics of the TT decay 
products. 


5 Summary and Conclusions 

The existence of tiny neutrino masses has broad impact in cosmology, particle, and nuclear 
physics. Hence, understanding the origin of their sub-eV masses is a pressing issue. We 
report the leading QCD corrections to the production rate and kinematic distributions of 
hypothetical, heavy Type III Seesaw lepton pairs, including soft gluon resummation of the 
heavy dilepton system. We find: 

1. The pp —>• and T+T“ iC-factors at NLO in QCD (see table 1) span 

: 1.17 - 1.37 (1.17 - 1.36) [1.14 - 1.17] at ^/s = 13 (14) [100] TeV(5.1) 
r+T" : 1.18 - 1.35 (1.19 - 1.33) [1.15 - 1.29] at ^/s = 13 (14) [100] TeV(5.2) 
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2. For the range of rriT considered, the total and T~^T~ production cross sections 

exhibit a scale dependence at 14 TeV and ±1% dependence at 100 TeV; see 
section 4.2. 

3. Differential iF-factors with respect to rapidity, transverse momentum, and polar dis¬ 
tributions of heavy triplet leptons are largely flat, indicating that naiVe rescaling of 
Born-level result is a largely justified estimate of kinematics at 0{as)', see section 4.3. 

4. The resummed transverse momentum of the dilepton system illustrates that 

the impact of TeV-scale systems recoiling off soft radiation is large; see figure 7. 

5. TT production and decay into the W + 4j final-state has been estimated at the HL- 
LHC, showing maximally a 4.8 — 6.3 cj discovery potential for rriT = 1.5 — 1.6 TeV. 
At 14 TeV and after 300 (3000) fb“^, there is 2cr sensitivity up to rriT = 1-3 — 
1.4 TeV (1.7— 1.8 TeV) in the individual CC and NC channels. At 100 TeV and after 
10 fb“^, a 5cJ discovery can be achieved for rriT ~ 1-4 — 1.6 TeV. See section 4.4. 
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A Triplet Lepton Pair Production at NLO in QCD via PSS 

The PSS technique and explicit examples, including the DY process, are reviewed authori¬ 
tatively in Ref. [36] . Terms relevant to triplet lepton pair production in hadron collisions are 
collected here for convenience. We emphasize that the NLO calculation greatly simplifies 
to simply convolving PDFs with tree-level amplitudes. 

A.l QCD-Corrected Two-Body Final State 

In the soft/collinear limits, radiation in TTj production becomes unresolvable and the 
three-body kinematics approach those of the two-body process 

(1{pa) q'ijPB) -t T{pi) T{p2). (A.l) 

Subsequently, the singular propagators in the amplitude factor and the soft/collinear contri¬ 
butions to the TTj cross section can be expressed as a divergent, but process-independent, 
piece and the LO TT cross section. Loop corrections cancel soft and soft-collinear poles 
rendering finite the quantity 

d(T{2) = d&^ + da^ + da^ + da^^. (A.2) 

PDFs are redefined to subtract residual hard-collinear divergence. We now give each term 
explicitly. 
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A.1.1 Virtual Corrections 


For massless fermions, 0{as) corrections to EW vertices, i.e., figure 1 (b), factorize [61] 
into a product of the Born matrix element, Mb-, and a universal form factor: 

Hpd)l^ {olPl + gnPR) u{pu) ^ HPd)'i^^{Pu,Pd)u{pu), (A.3) 


v{pd)T^{Pu,Pd)u{pu) = v{pd)'y^ [qlPl + OrPr) u{Pu) x F 
_ as{pl) ^ 


F = 


CM = 


-CirC,(s)(-irr(l + e)r(l-e) 


dvr 

r (1 — e) 


-2 


r (1 - 2e) ^ 


Cf = 4/3. 


(A.4) 

,(A.5) 

(A.6) 


The vertex correction is UV-finite, so its counter term is zero. Counter terms cancel external 
quark self-energy corrections in the on-shell and MS schemes. Hence, the summed, squared 
amplitude is 


[TWIjTW i_Loop] + 0{al), 7Wi_Loop = MbF (A.7) 
= J^I-^bP(1 + 25R[.T])+C)(a2). (A.8) 


The Born and virtual terms at 0{as) for heavy lepton production are then 

- 2 ^ AY 


da^ -I- dd'' = dd 


V 




1 + ^4^CfCs(s) 


2f 


fAY AY 




AM -2, A=-3, Ao"=-8 + —. 


(A.9) 

(A.IO) 


A.1.2 Soft Corrections 

For colored partons a,b, and color-connected Born cross section the soft radiation 
expression is 


s _ OIs{pY^ 


dd^ = 


2f 


Ceis) [ 

^ U — ^ 7T' 


dS 


a^b=q,q' 


-jPa-Pb) 

{Pa- Pj){Pa- Pj)' 


= {-i)^4^cm [ ds 


2f 


{Pq-Pg') ^.0 ^ jPg' ' Pq) 


ddC, 


dS = 


1 /4 


F \S 


L 


SsVS/2 


dEjE]-^ 


.{Pq-PjY {Pq'-PjY '''' 

d9suY~‘^'^6 J d(psm~‘^^ 4>, 


(ATI) 

(A.12) 

(A.13) 


where dS is the soft one-particle phase space in n = 4 — 2e dimensions, and we make use 
of the fact that the only two colored partons in the DY process are the initial-state quark 
and antiquark. For the present case, the color-connected and plain Born cross sections are 
related by 

ddj-, = ddjg = -Cf dd^. (A. 14) 

We take Pg (pgi) to propagate in the Fz (—z) direction, implying 


VI 


VI- 


Pq ■ Pj = ^Pj(^ - cos 9), pg^-pj =—Ej(1Fcos 9), pg-p^, = -. (A.15) 
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da^ = da^ 

\S 


Using the tabled integrals of Ref. [36], the soft contribution to DY processes is 

-^CpC.is) ^_ + _ + AojJ, 

^2 =2, = -dlogJs, ^0 = 41og^(5s. 

A.1.3 Soft Collinear Corrections 

For qq' initial-states, the soft-collinear correction to the Born amplitude is 
asifJ-r). 


(A.16) 

(A.17) 


da^^ = da 


B 


2-k 


-CpCeis) 


X 


= da 




qg) 


2-k 


CpCeis) ( —Y- - + 2AQ^(g —)• qg) 


Af^ = ( 2 log ^5 + 0 , = 1^2 log ^5 + 0 log 


A.18) 

(A.19) 

(A.20) 


The soft-collinear splitting functions, A(1 —)• 2 3), are derived by integrating over the 
Altarelli-Parisi (AP) splitting functions. The <?—)■(? and q' —>• q' functions are equal by 
CP-invariance and therefore are combined in the last step. 

A.1.4 Assembly of Two-Body Corrections 

A finite result emerges after summing the virtual, soft, and soft-collinear corrections. 


da^ -|- da^ + dd^^ = da 


sc 




27r 


-CpCeis) 


A 2 , ^ 2 ^ , ^ Af ^ 2Af^(q^qg) 

e e 


2 “I— Y I 
e"' / V e 


+ (A^ + Aq -|- 2AQ^(q —)• qg)) 


(A.21) 


and using Eqs. (A. 10), (A.17), and (A.20) to show 

-b Af = 0, AX + Af + Af^ = 0. (A.22) 


With poles removed, we send e ^ 0. The non-hard-collinear, two-body contribution of 
Eq. (3.10) is 


'.B 


d(T(2) = + dd^ + dd° + dd‘ 


V 




= dd 


B 




(A.23) 
(A.24) 


A.2 Hard-Collinear Subtraction 

Remaining hard-collinear ISR poles are associated with PDFs, which themselves are all or¬ 
ders resummations of hard-collinear splittings. The divergences are removed by subtracting 
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out from PDFs redundant 0{as) splittings. The hard-collinear subtraction is given by 


a 


HC/ 


\pp^TTX)= Y, [ f di2 &^{ab ^TT) 

a,b=q,q "^o/Cl 

X /a/p(6,/^^)/b/pfe,/^^) +/a/p(6,/^^)/b/p(6,/^^) + (1 ^ 2|.\,.25) 
where the redefined, 0(as)-corrected PDFs /a/p(C)/^^) ^-^e 

(A.26) 




27r 


b=g,q,g' ^ 


fb/p 1 11 Pab\^) 


j 2-^g/p + y ^/p/p 


The summation is over all partons that give rise to q, q' after one radiation, 6ba is the 
Kronecker d-function, and the modified AP i ^ j splitting functions Pji{z) have the form 


Pji{z) = Pji{z) log 


(1 - z) 5cs 


-P'niz). 


z 

In n = 4 — 2e dimensions, P and P' are related by Pji(z,£) = Pji{z) + ePj^{z), with 

(1 + ^2) 


(A.28) 


Pqqi^) = Cp- 


1 - Z 


P', = -C'i.(l-Z), 


1 


Pqgi^) = 9 + (1 - ’ Pqg(^) = -^(1 - 


(A.29) 


A.3 Efficient Event Generation of Three-Body, Hard, Non-Collinear Radiation 

Imposing soft/collinear cuts in the partonic c.m. frame on the processes 

q{paW{pb), g{PA)q{jPB), g{PA)q'{pB) -t T{pi) T{p2) j(pj) (A.30) 

regulates all diverges. TTj event generation can be handled efficiently by implementing 
Eqs. (3.3)-(3.5) into MGS. This entails inserting into the header of the file SubProcesses/cuts. f: 

double precision del_s, del_c, sHat, ejMin, tjMin, tmptj 

and in the file’s body: 

del_s = 3.0d-3 

del_c = del_s * l.Od-2 

sHat = SumDot(p(0,1),p(0,2),IdO) 

ejMin = del_s * dSqrt( sHat ) /2.0d0 

tjMin = del_c * sHat 


do i=l.nincoming 

do j =nincoming+l,nexternal 
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if(is_a_j(j)) then 


if( p(0,j).LT.ejMin ) then 

passcuts=.false. 

endif 

tmptj = dabs(SumDot(p(0,i),p(0,j),-ldO)) 
if( tmptj.LT.tjMin ) then 
passcuts=.false. 
endif 


endif 

enddo 

enddo 

For TTj, Sij cuts are not needed as all s-channel poles are regulated by nonzero rriT- 


B Recoil qt Resummation with Fixed Order Matching 


Individual formulae and steps for resumming recoil logarithms at LL accuracy for DY-type 
processes within the CSS formalism [37-39] can be found across literature [39, 51, 52, 62-67] 
For conciseness and convenience, they are collected here in terms of generic tree-level 
quantities. 

For the processes 

qq'^V*^TT, Ye{7,Z,VF}, (B.l) 

the resummed and matched qt distribution of the TT system decomposes into the FO, 
resummed, and asymptotic pieces; see Eq. (3.17). The oscillatory behavior of Jo{qTb) and 
the fine numerical cancellation of the resummed and asymptotic terms when qt > My* give 
rise to an unnecessary computational burden for a contribution that is small compared to the 
FO piece. In practice [62, 63], one introduces that is unity in the qx/My* —)■ 0 

limit and vanishes the asymptotic-resummed difference at qx > My*\ 

l + (B-2) 


We set g^B,tch _ y[y, [ 52 ]^ and write the matched, triply differential distribution [51, 52] 


da 


Matched 


da™ 


dMyt dy dq^ dMy, dy dqi 


+ f 


Match 


(qt) 


^^Resum 


Cio-Asymp 


dMyt dy dq^ dMy* dy dqj 


. (B.3) 


In the qx/My* —?■ 0 limit, the longitudinal momentum fractions of initial-state partons q, q' 
are 


ei = e^ 


My* 


^2 = e y 


My* 

7^’ 


dMy* dy 


s dii di2- 


(B.4) 


^We report and clarify a slight but unfortunate typo in Ref. [52]: the resummed result is to LL accuracy, 
not NLL. 
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Though untrue for the resummed and asymptotic pieces at large qx/My*, their combined 
contribution vanishes by construction. For the FO calculation in this limit, My* —?> \/I and 
the above momentum fractions are for initial-state partons q,q', and g. With this, one may 
write 


d^A ^ ^^Resum 

dMyt, dy dq^ dMy, dy dq^ dMy^^ dy dq^ 


1 / (i^Resum 

s \d^i d ^2 dq^ 


dfjAsymp 

d^2 dq^ 


(B.5) 


Dehning da^/d^id^ 2 dqT as the quantity in the parentheses, one at last has 
^^Matched ^ rl rl da^ _ 4m| 

dqr dqx Xo ^ Xq/^i ^ <^6 d ^2 dqr ’ s 


(B. 6 ) 


To account for the appropriate normalization at 0{as), the area under the above distribu¬ 
tion is scaled to equal in Eq. (3.11). The resummed and asymptotic terms are now 

discussed. 


B.l Recoil qx Resummation 


Recalling the resummed expression, but in terms of 

j_Resum -i roo 

,, ,, , 2 = 7 / db^ MbQT) X W X a^iqq^T T), (B.7) 

d^i d^2 dq^ 4 

W = ^-Sp{b„Mv*) W(6*,6,6). (B.8) 


Technically, W is ill-defined for b > 1/Aqcd due to divergent as{b~‘^). This is amelio¬ 
rated [39, 64] by Sjyp, which is valid for all b and approaches unity when b <C I/AqcDj and 
by introducing 


b,{b) 


b 

\/l + ^^/^max 


< b 


max? 


1/6 

max > Aqcd- 


(B.9) 


This approach, however, requires that Sjyp be extracted from data and dependent on 6 inax- 
We choose the non-perturbative Sudakov factor given by the BLNY parameteriza¬ 
tion [65, 66 ] 


SNp{b,Q,^ 1 ,^ 2 ) 


91 + 92 log 


Q 


2Qblny 


+ 5'i53log(100^i6) 


(B.IO) 


Fits for Qi to Tevatron Z boson data with Qblny = 1-6 GeV and 6max = 0.5 GeV ^ 
give [ 66 ] 

gi = 0.211°;°} GeVl 52 = 0.681°;°} GeVl 53 = -0.60l°;°^ (B.ll) 

To numerically integrate over the domain of 6 , we exploit the 6 -dependence oi Snp and 
dehne: 


6,(6) = e 


= 


Snp — SNP—b'^9i = 


92 log 


Q 


2(5blny 


+ 5153 log(lOO^ACB) 


. (B.12) 


The 6 integral now takes the manageable form 



g-Sjvp 


1 

51 


dh 


(B.13) 
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The perturbative Sudakov-like form factor is given by 


Jcl/hl 


Alog^ + B 


(B.14) 


where the low-scale integration limit is /r = Co/6*, with cq = 2e~^^ and '^e ~ 0.577 is 
the Euler-Mascheroni constant. The functions A and B can be expanded perturbatively in 
powers of For a resummed calculation at LL, we take A and B to The 

expressions are [39] 


•4 = E-4'”>(£) =-4'‘>(S)+0(a?), A<‘) = 2Cp, 


n=l 

oo 


2tt 


® = E (S)" = (S) + = -3Cf. 

n=l 

At one-loop, Og is given by 

= - 




1 +bo as{Ml) log 


ill 

Ml 




and allows Eq. (B.14) to be evaluated analytically [67]. The result is 


Sp{b^,Q) 


rQ2 


dl OLsil^) 


cUbl d- 

Cf 




0 1 <5" 

2 log ^ - 3 


2-Kb^a]{M^) + 26 oa.(M|)t(Q 2 ) _ 36o«.(M|)] log 

'blQl 


X 


Ce^ 

— 

vr&o 


A = 


1 boasiM^)t{Q^) 

1 6oa^(M|)t(cg/62)' 


V) = log^. 


(B.15) 

(B.16) 


(B.17) 


(B.18) 


(B.19) 

(B.20) 


In the perturbative limit b <C I/AqcDj the TMD factorize into universal Wilson co¬ 
efficient functions for i ^ j splitting, Cji, and the usual transverse momentum-independent 
PDFs: 


^q/piCl,b^, lf)J^lp{C2,b‘^, If) = ^ [Cqi'Sl fi/p] {Cl, If) X [Cq'k'^ fk/p] (6,/^/)- 

i,k=q,q',g 

(B.21) 

The summation is over all partons i,k that can split into q,q'■ The convolution notation 
denotes ^ 

[Cji ® h/p] iC, ^ y C,i{z) h/p 0, (B.22) 

Like A, B in Sp, Cji can be expanded in powers of as and we take Cji to 0{a'^^^) [39]: 

CX) 

C,,{z) = Y,C^\z) « C-]?(z) = 6,^1 - z). (B.23) 

n=0 
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Formally, Cji{zY^^ is an 0{ag) contribution except at qt = 0, where it 0(ag). Hence, our 
qt spectrum is accurate to 0(as) in shape everywhere but the origin [52], an acceptable 
error considering the uncertainty in our knowledge of rriT- Evaluating the convolutions, the 
product of TMDs is 


~ Cj°'> 0 fq/p X ® fq/p = /g/p(6, /^/)/g/p(6, M/)- 


Equating the collinear and impact scales, = Cq/6 ^, 1F(6 *,^i,^2 ) becomes 


(B.24) 
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fq/p ( Cl) ^2 ] fl/P [ CS) q,2 ] fl/P [ ^2) ^2 ] fl/P [ 


hi 


hi 


(B.25) 


Though valid only for qx/My* —)• 0, extension to 6 > Aqcd is allowed hy Sjyp and Jo{x), 
which vanish for large b and x = qxh. After simplification, the resummed expression to LL 
accuracy is 


^^Resum 


= ^ I dh Mbqx) e-Sp{b„Mv*) W{b^,^ 1 ,^ 2 ) x a 


dCi d^2 dqx 2gi Jq 


B 


(B.26) 

where expressions for the relevant factors are given in Eqs. (B.12), (B.19), and (B.25). In 
essence, the expression preceding is the unintegrated TMD parton luminosity for DY 
systems. 


B.2 Asymptotic Expansion of Resummed Expression 

The asymptotic piece is obtained by formally expanding the resummed expression in powers 
of Og, keeping only terms as singular as l/q^, and evaluating the impact parameter integral. 
The result is 


^^Asymp ^ ^^^^2) 


E 


A\0g^^ + B \ fq/p{£.l,ljf)fq/p{^2,d^) 
Qt / 


dCi d ^2 dqx qr 'x 

yiy L 

+ fq/pi^l^lf^)fq/p{^2,lf^) + fq/p{^l, d‘^)fq/p{C2, ff^) + (Cl ^ C 2 ) 


X d^B.27) 


The renormalization and factorization scales here must be set equal to the FO scales, i.e., 

Eq. (4.2), in order to avoid spurious logarithms containing ratios of factorization/renormalization 
scales. To 0{as), A and B are given in Eqs. (B.15)-(B.16). The PDFs, corrected for single 
parton splitting, are 


/g/p(C,fi^) 


2-k 

2-k 


[Pgi^fi/p] (C,A^) 

^=1,9 

f T A/f 

Ui ^ 


(B.28) 
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where the AP splitting functions are in Eq. (A.29) and the “plus distribution” [P(z)]_|_ is 
dehned as 


f 


dz [Piz)]+ /( 


dz [P{z)]+ /( 



dz P{z) [f{z,fi‘^) - f{l)], and (B.30) 

dz P{z) [/(z,/i^)-/(!)] +/(I) / dz P{z). (B.31) 

Jo 


Expanding the q ^ q integral under the plus distribution gives the expression 


J ^ fq/p P 


and allows one to write 




1 




d^ ^^qqi^) fq/p d' J _ ^'jfq/pi^^ d ) 


2Cf 


+ 2 [1 + Ci? log(l — ^)] /q/p(^,/^^), 


(B.32) 


fq/piC, d)) = ^ y Pqq{z) fq/p d^) + Pqg{z)fg/p d' 


+ 


2tt 
asilJ?) 


vr 


[1 + Cf log(l - 0] /g/p(?: d^)^ 


(B.33) 


There is slight technical distinction between / here and / in Eq. (A.27), where hard collinear 
splittings in PSS are addressed. The / are regulated by subtracting out individual singular 
points via the plus distribution, whereas / are regulated by the soft cutoff 6s- 
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